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^ i Abstract 

,—1. \ Based on a scenario of the inhomogeneous big-bang nucleosynthesis (IBBN), we investigate the 



detailed nucleosynthesis that includes the production of heavy elements beyond ^Li. From the 
observational constraints on light elements of ^He and D for the baryon-to-photon ratio given by 
WMAP, possible regions found on the plane of the volume fraction of the high density region 
against the ratio between high- and low-density regions. 



X 

C^ ■ In these allowed regions, we have confirmed that the heavy elements beyond Fe can be produced 

appreciably, where p- and/or r-process elements are produced well simultaneously compared to the 
solar system abundances. We suggest that recent observational signals such as ^He overabundance 
in globular clusters and high metallicity abundances in quasars could be partly due to the results 
of IBBN. Possible implications are given for the formation of the first generation stars. 

PACS numbers: 26.35.+C, 98.80.Ft, 13.60.Rj 
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I. INTRODUCTION 



Big bang nucleosynthesis has been investigated mainly on the context of the standard 
cosmological model (SBBN), where origin of light elements of ^He, D, and ^Li have been 
discussed in detail [l|. While observations of ^He are still in debate with the uncertainty 
of 20-30 % in the abundance (2t^, those of D constrain severely the possible range of the 
abundance production in the early universe p-Ll!]. Contrary to the above standard BBN, 



the heavy element nuc 
from twenty years ago 



eosynthesis beyond the mass number A = 8 has been proposed 
8l-ll4|. where the model is called the inhomogeneous BBN (IBBN). 
This model relays on the inhomogeneity of baryon concentrations that could be induced 
by baryogenesis (e_.g; ^Ref. IJ]) or some phase transitions such as QCD or electro- weak 



phase transition [8 
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16| during the expansion o 
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the universe. Although a large scale 



181], small scale one has been advocated 



inhomogeneity is inhibited by many observations 

within the present accuracy of the observations. Therefore, it remains a possibility for IBBN 

to occur in some degree during the early era. 

On the other hand, Wilkinson Microwave Anisotropy Probe (WMAP) has derived crit- 
ical parameters concerning the cosmology of which the present baryon-to-photon ratio rj is 
determined to he t] = (6.19 ± 0.15) x 10^^° 18|. This value is almost consistent with that 
obtained from the observation of D. Therefore, considering the uncertainty of the ^He abun- 
dance, we can fix the ratio t] in the discussion of the nucleosynthesis in the early universe. 
If the present ratio of t] is determined, BBN can be performed along that line in the ther- 
modynamical history with use of the nuclear reaction network. On the other hand, peculiar 
observations of abundances for heavy elements and/or ^He could be understood in the way 
of IBBN. For example, the quasar metallicity of C, N, and Si could have been explained 
from IBBN [19]. Furthermore, from recent observations of globular clusters, possibility of 



inhomogeneous helium distribution is pointed out 



where some separate groups of dif- 



ferent main sequences in blue band of low mass stars are assumed due to high primordial 
helium abundances compared to the standard value 20|, |21 |. 

Despite a negative opinion against IBBN due to insufficient consideration of the scale 
of the inhomogeneity 22J, Matsuura et al. have found that the heavy element synthesis 
for both p- and r-processes is possible if r^ > 10~^ 23|, where they have also shown that 
the high rj regions are compatible with the observations of the light elements, ^He and 



D [2J]. However, their analysis is only limited to a parameter of a specific baryon number 
concentration. Therefore, it should be needed to constrain the possible regions from available 
observations in the wide parameter space that describes the IBBN. 



In ^mi we review and give the adopted model of IBBN 2J]. Constraints on the critical 
parameters of IBBN due to light element observations are shown in §111, and the productions 
of possible heavy element nucleosynthesis is presented in §IV. Finally, ||V] is devoted to the 
summary and discussion. 

II. COSMOLOGICAL MODEL 

We adopt the two-zone model for the inhomogeneous BBN, where the early universe 
is assumed to have the high- and low- baryon density regions [12] under the background 
temperature T. For simplicity we ignore the diffusion effects before (10^°K < T < lO^^K) 
and during the primordial nucleosynthesis (lO^K < T < 10^'^K), where the plausibility will 
be discussed in ^|Vl After the epoch of BBN, all the elements are assumed to be mixed 
homogeneously. 

Let us define the notations, n^^^,nf^^g,^, and n^^^ as averaged-, high-, and low- baryon 
number densities. /^ is the volume fraction of the high baryon density region. Xf*"^,Xj *^ 
and Xl°^ are mass fractions of each element i in averaged-, high- and low-density regions, 
respectively. Then, basic relations are written as follows: 






^ave^i — tv'^hiqh^i + (^ tv) "^low-^i ■ (2) 



Here we assume the baryon fluctuation to be isothermal as was done in previous studies (e.g.. 



M 



Refs. [8|, y, ll5|). Under that assumption, since the baryon-to-photon ratio is defined by the 
number density of photon in standard BBN, n^ = 2^(3)/7r^ [ksT /he) , Eqs.([T]) and ([2]) are 
rewritten as follows: 

Vave = fvVhigh + (^ ~ fvjVlow^ (3) 

high 



Va.eXr = LX^'Sm^H + (1 - QX^'^Vlo^. (4) 



present 



where r^s with subscripts are the baryon-to-photon ratios in each region. In the 
paper, we fix rj^^^ = 6.1 x 10^^° from the cosmic microwave background observation |l7l. Il8| 



rjf^igf^ and 77^^^ are obtained from both /^ and the density ratio between high- and low-density 
region: R = n^igh/n^^^ = r^ughhiow 

To calculate the evolution of the universe, we solve the following Friedmann equation, 



xV SttG 



P. (5) 



^x J 3 

where x is the cosmic scale factor and G is the gravitational constant. The total energy 
density p is the sum of decomposed parts: 

P = Py + Pe±+ Pu + Pb- 

Here the subscripts 7, e"^, u, and b indicate photons, electrons/positrons, neutrinos, and 
baryons, respectively. We note that p^ is the average value of baryon density obtained from 
Eq. dH). 

The energy conservation law is used to get the time evolution of the temperature and the 
baryon density, 

|(,xV4(.»)^0. (6) 

where p is the pressure of the fluid. 

III. CONSTRAINTS FROM LIGHT-ELEMENT OBSERVATIONS 

In this section, we calculate the nucleosynthesis in high- and low-density regions with use 



of the BBN code 27 which includes 24 nuclei from neutron to ^^O. We adopt the reaction 
rates of NACRE 



311], the neutron life time r^ = 885.7 sec J33i], and take account of the 
number of species of the massless neutrinos N^^ = 3. 

Figure [1] illustrates the light element synthesis in the high- and low-density regions with 
/^ = 10-6 and R = 10^ that correspond to r/^^^^ = 3.05 x 10"^ and r/,„^ = 3.05 x 10-^°. In 
the low-density region the evolution of the elements is almost the same as that of standard 
BBN. In the high-density region, while ^He is more abundant than that in the low-density 
region, ^Li (or ^Be) is much less produced. It implies that heavier nuclei such as ^^O, hardly 
synthesized in SBBN, are synthesized at high-density region. 

For /^ <^ 0.1, the heavier elements can be synthesized in the high-density regions as 



discussed in Ref. 



13|. For /„ > 0.1, contribution of the low-density region to rj^^^ can be 



neglected and therefore to be consistent with observations of light elements, we need to 
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Fig. 1: Illustration of the nucleosynthesis in the two-zone IBBN model with /^ = 10~^ and R = 10^. 
The baryon-to-photon ratios in the high- (left panel) and low- (right panel) density regions are 
''Ihigh ~ ^-^^ ^ ^^'^ ^^*^ ''llow ~ ^-^^ ^ 10^^°, respectively. 



impose the condition of /^ < 0.1. Now, we put constraints on /^ and R by comparing the 
average values of ^He and D obtained from Eq. (j4]) with the following observational values. 
First we adopt the primordial ^He abundance reported in Ref. 



4: 



0.232 <Yp< 0.258. 



(7) 



Next, we take the primordial abundance from the D/H observation reported in Ref. 28 1 



D/H= (2.84 ±0.26) X 10-^ (8) 

where the systematic error given in Ref. [7| is adopted. 

Figure [2] illustrates the constraints on the f^ — R plane from the above light-element 
observations with contours of constant rj^igh- The solid and dashed lines indicate the upper 
limits from Eqs. ([7]) and ([8]), respectively. As the results, we can obtain approximately the 



following relations between f^ and R : 

0.26 X /-0-96 for /^ > 3.2 x 10-^ 
R<{ (9) 

1.20 X /-°-83 fQ^ j^ < 3 2 X 10-^ 

As shown in Figure [21 we can find the allowed regions which include the very high-density 
region such as %jg^ = 10~^. 

Matsuura et al. \2^ defined a parameter of the baryon number concentration a in the 
high density region instead of two parameters of /„ and R that are needed to solve Eqs.([3]) 
and dH): 

fvVhigh ■ (1 - fv) Viow = a : (1 - a) . 

However, they have only examined the case of %jg;i = 10~^ and 77^^^ = 3.162 x 10^^°, where 
a = 0.48 for rjave = 6.1 x 10^^°. Our constraints in Eq. (^ correspond to a = 0.02 — 0.5. 
Since we have fixed the value of rjave, we can obtain physically more reasonable regions on 
the plane of {f^, R). 

Naturally, as rjf^^, takes larger value, nuclei which are heavier than ^Li are synthesized 
more and more. Then we can estimate the amount of total CNO elements in the allowed 
region. Figure [3] illustrates the contours of the summation of the average values of the 
heavier nuclei {A> 7), which correspond to Fig. [2] and are drawn using the constraint from 
'^He and D/H observations . As a consequence, we get the upper limit of total mass fractions 
for heavier nuclei as follows: X{A > 7) < 10^^. 

We should note that abundance flows proceed beyond the CNO elements thanks to the 
larger network for high r/- values as shown in Table [TTl of the following section. 

IV. HEAVY ELEMENT PRODUCTION 

In the previous section, we have obtained the amount of CNO elements produced in the 
two- zone IBBN model. However, it is not enough to examine the nuclear production beyond 
A > 8 because the baryon density in the hi gh- d ensity region becomes so high that elements 



beyond CNO isotopes can be produced [lO|, 
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34| 



In this section, we investigate the heavy element nucleosynthesis in the high-density region 
considering the constraints shown in Fig. [2l The temperature and density evolutions are the 
same as used in the previous section. Abundance change is calculated with a large nuclear 
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Fig. 2: Constraints on the fy — R plane from the observations of hght element abundances. The 
region below the solid line is allowed one obtained from ^He observation 4|. Constraints from the 
D/H observation [28] are shown by the region below the dot-dashed line. The shaded region is the 
allowed parameters determined from the two observations of '^He and D/H. The dotted lines show 
the contours of the baryon-to-photon ratio in the high-density region. Filled squares indicate the 
parameters for heavy element nucleosynthesis adopted in ^IV[ 



reaction network, which includes 4463 nuclei from neutron {n), proton {p) to Americium 
{Z = 95 and A = 292). Nuclear data, such as reaction rates, nuclear masses, and partition 
functions, are the same as used in [26| except for the neutron-proton interaction; We take 
the weak interaction rates between n and p from Kawano code (29|, which is adequate for 
the high temperature epoch of T > 10^° K. We note that mass fraction of ^He and D 
obtained with the large network are consistent with those in in §1111 within the accuracy of 
few per cents. 

As seen in Fig. |3l heavy elements oi X{A > 7) > 10^^ are produced nearly along the 
upper limit of R. Therefore, to examine the efficiency of the heavy element production, we 
select five models with the following parameters: rjj^^ , = 10"'^, 5.3 x 10""^, 10""^, 5.5 x 10^^ 
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Fig. 3: Contours of the averaged total mass fractions which are the sum of nuclei heavier than ^Li, 
where we find consistent regions with ^He and D observations. 



, and 10-5 corresponded to {f^,R) = (2.1 x 10-^ 1.8 x 10*5) , (7.0 x 10-^ 9.3 x 10^) , 
(5.7 X 10-^1.9 X 10^) , (1.0 X 10-^9.6 x 10^) , and (7.5 x 10-^,1.9 x 10^). Adopted pa- 
rameters are indicated by filled squares in Fig. [2l 

Figure H] shows the results of nucleosynthesis in the high-density regions with rjf^-f^ ~ 
10"^ and 10"^. For rjfj^-f^ ~ 10"^, the nucleosynthesis paths are classified with the mass 
number [23|. For nuclei of mass number A < 100, proton captures are very active compared 
to the neutron capture of T > 2 x 10^ K and the path moves to the proton rich side, which 
began by breaking out of the hot CNO cycle. For nuclei of 100 < A < 120, the path goes 
across the stable nuclei from proton to neutron rich side, since the temperature decreases and 
the number of seed nuclei of the neutron capture process increase significantly. Concerning 
heavier nuclei of ^4 > 120, neutron captures become much more efficient. In Figure ll](a), 
we see the time evolution of the abundances of Gd and Eu for the mass number 159. First 
i59rpj^ (stable r-element) is synthesized and later ^^^Gd and ^^^Eu are synthesized through 
the neutron captures. After t = 10^ sec, ^^^Eu decays to nuclei by way of ^^^Eu — > ^^^Gd 
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_s, i59rp|^^ where the hfetimes of ^^^Eu and ^^^Gd are 10.1 min and 18.479 h, respectively. 
These neutron capture process is not similar to the canonical r— process, since the nuclear 
processes proceed under the condition of the high-abundance of protons. 

For r]f^- f^ ~ 10^"^, the reactions first proceed along the stable line, because triple-a re- 
actions and other particle induced reactions are very effective. Subsequently, the reactions 
directly proceeds to the proton rich region, through rapid proton captures. As shown in 
Fig. m^b), ^''^Sn which is proton-rich nuclei is synthesized. After that, stable nuclei ^°^Cd 
is synthesized by way of ^°^Sn -^ ^°^In -^ ^'^^Cd, where the lifetimes of ^°^Sn and ^"^In are 
10.3 min and 58.0 min, respectively. In addition, we notice the production of radioactive 
nuclei of ^^Ni and ^'^Co, where ^^Ni is produced at early times, just after the formation of 
■^He. Usually, nuclei such as ^^Ni and ^^Co are produced in supernova explosions, which are 



assumed to be the events after the first star formation (e.g. Ref. 32|). In IBBN model, how- 
ever, this production can be found to occur at extremely high density region of rj^^-gf^ > 10^^ 
as the primary elements without supernova events in the early universe. 

To explain differences of the nuclear reactions which depend on the baryon density, we 
focus on the neutron abundances. Figure [5] shows the evolutions of the neutron abundances 
in the SBBN and IBBN models. For rjhigh = 10~^, neutron abundance decreases rapidly 
at ~ 10 sec to the formation of ^He and ^^Ni. Thus, neutron abundance is not enough to 
induce the neutron capture producing heavy nuclei of A > 90. On the other hand, neutron 
abundance tends to remain even at the high temperature for the lower value of ?7/ij„/j. We 
can see the case of f]hi„h = 10~^, where there remain much neutrons to occur the neutron 
capture reaction. Thus the neutron capture process to produce heavy elements of A > 90 
can become active. 

Time scales in the decrease for the neutron abundances change drastically the flow of the 
abundance production. Figures |6] and [7] show the flows for r^^j^^ = 5.3 x lO"''. Before the 
signiflcant decrease in the neutron abundances before 10 sec, the nucleosynthesis proceeds 
already along the stable line by way of the neutron included reactions (Fig. |6]). At that 
time, the nuclear reactions are stuck around Z = 60 with A^ = 82, since it takes time to 
synthesize heavier nuclei because Nd {Z = 60) and Sm {Z = 62) have some stable isotopes. 
As time goes, neutron captures of these nuclei start, where the neutron captures proceed 
signiflcantly and r-elements can be synthesized. After the depletion of neutrons {t > 40 sec), 
nuclei around the neutron numbers A^ = 82 are produced through proton induced reactions 



such as ^^^Sm (Fig. [7]). 

Final results (T = 4 x 10^ K) of nucleosynthesis calculations are shown in Tables HI and HTl 
Table H] shows the abundances of light elements, ^He, D, and ^Li, in high- and low-density 
regions with their average values. Abundances of the low-density side (the third and sixth 
columns) are obtained from the calculation by BBN code used in §111, because abundance 
flows beyond A = 7 are negligible. We should note that the average abundances of ^He and 
D are consistent with their observational values of ([7]) and (|8]) . Table [TTl shows the amounts 
of heavy elements. When we have calculate the average values, we set the abundances of 
y4 > 16 as zero for low-density side. For rji^- , ~ 10~^, a lot of nuclei oi A > 7 are synthesized 
whose amounts are comparable to that of ^Li. Produced elements in this case include both s- 
element (i.e. ^^^Ba) and r-elements (for instance, ^^^Ce and ^^^Nd), since moderate amounts 
of neutrons remain as shown in Fig. |5] 

For rjf^-f^ ~ 10~^, there are few r-elements while both s-elements (i.e. ^^Kr and ^^Y) 
and p-elements (i.e '^^Se and '^^Kr) are synthesized such as the case of supernova explosions. 
Although heavy nuclei oi A > 100 are not synthesized appreciably, those oi A < 90 are 
produced well owing to the explosive nucleosynthesis under the high density circumstances 
(p ~ 10^ g cm~'^). The most abundant element is found to be ^^Ni whose production value 
is much larger than the estimated upper limit of the total mass fraction (shown in Figl3]) 
derived from the BBN code calculations. This is because our BBN code used in §1111 includes 
the elements up to A = 16 and the actual abundance flow proceeds to much heavier elements. 

Figure E] shows the abundances averaged between high- and low-density region using 



Eq. (jlj) compared with the solar system abundances [30|. For rjf^^gf^ ~ 10~^, abundance 
productions of 120 < A < 180 are comparable to the solar values. For rii^-^f^ ~ 10"'^, those 
of 50 < A < 100 have been synthesized well. In the case of ^^^^/j = 5.3 x 10~^, there are 
outstanding two peaks; one is around A = 56 {N = 28) and the other can be found around 
A = 140. Abundance patterns are very different from that of the solar system, because 
IBBN occurs under the condition of significant abundances of both neutrons and protons. 

V. SUMMARY AND DISCUSSION 

We have investigated the consistency between inhomogeneous BBN and the observation 
of '^He and D/H abundances under the standard cosmological model having rj determined 
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Table. I: Mass fractions of light elements for follow cases : r/^. 



igh 



w 



Vhigh 



5 X lo~^ 



r/^. ^ ~ 10 ■*, and f?;^j„/j = 10 ^. t r^^ and Tr-^ is the time and temperature at the final stage of the 



calculations. 








(''Ihigh'Vlow) 


2.1 X 10-^1.8x 10^ 
(1.06 X 10-3, 5.88 X 10-10) 




7.0 X 10-^9.3x 10^ 
(5.33 X 10-^ 5.73 x lO-^O) 


yfin-'^fin) 


1.0 X lO^sec, 4.2 x 10^ K 




1.1 X lO^sec, 4.9 X 10^ K 


elements 


high low average 


high low average 



p 

D 

^He 

^Li 



0.586 

1.76 

0.413 



0.753 



0.746 



0.598 



0.753 



X 10-21 4.48 



0.743 
10-5 4.32 X 10-5 4.14 X lO-^i 4.67 x 10-^ 4.38 x 10-^ 



0.247 0.253 



: X 10-21 4.( 
0.402 



0.247 



0.257 



1.63 X 10-13 1.79 X 10-^ 1.72 x 10-^ 3.43 x lO-i^ 1.70 x 10-^ 1.59 x 10-^ 





(a) For cases of r]^-^^ = 10" 


~ 


and 


^fc 


igh 


= 5 X 10-1 


vlhigh^''llow) 


1.0 X 10-^1.7x 10^ 
(1.02 X 10-^, 6.00 X 10-1°) 










7.5 X 10-^1.9x 10^ 
(1.02 X 10-5, 5.34 X 10-1°) 


^fin-'^fin) 


1.2 X lO^sec, 4.3 x lO"^ K 










1.2 X lO^sec, 4.5 x lO"^ K 


elements 


high low avera 


,ge 






hi 


gh low average 



p 

D 

^He 

7Li 



0.638 



0.752 



0.750 



0.670 0.753 0.743 

6.84 X 10-22 4.34 x 10-^ 4.27 x 10-^ 1.12 x 10-22 5.19 x 10-^ 4.55 x 10-^ 



0.362 



0.248 



0.249 



0.330 



..-.. ..™. ..-_ 0.257 

,70 X 10-9 6.73 X 10-^ 1. 
(b) For " ^""^ ' ^""^ 



7.42 X 10-13 1.87 X 10-9 1. 

cases of r/^^^^ = IQ-^ and r/^^^^^ = 10" 



0.246 
.47 X 10-9 9.63 X 10-9 



by WMAP. We have adopted the two- zone model, where the universe has the high- and low- 
baryon density regions at the BBN epoch. 

First, we have calculated the light element nucleosynthesis using the BBN code having 
24 nuclei for the high- and low-density regions. We have assumed that the diffusion effect 
is negligible. There are significant differences for the time evolution of the light element 
between the high- and low-density regions; In the high-density region, the nucleosynthesis 
begins faster and ^He is more abundant than that in the low density region as shown in 
Figure HI From 'iRe and D/H observations, we can put severe constraint on two parameters 
of the two-zone model: the volume fraction /^ of the high-density region and the density 
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ratio R between the two regions, where we have assumed that abundances in the two regions 
are mixed homogeneously. 

Second, using the allowed parameters constrained from the light element observations, we 
calculate the nucleosynthesis that includes 4463 nuclei in the high-density regions. Qualita- 
tively, results of nucleosynthesis are the same as those in Ref. [2^]. In the present results, we 
showed that p- and r-elements are synthesized simultaneously at high-density region with 
r]high — 10~'*.d Such a curious site of the nucleosynthesis have never been known in previous 
studies of nucleosynthesis. 

As the results, we have obtained the average values of mass fractions from the nucleosyn- 
thesis in high-density and that in low density regions. The total averaged mass fractions 
beyond the light elements X{A > 7) are constrained to be 10~^ (for ^j^^j = 10"'^) and 
10~^ (for rjf^-^ = 10~^). We find that the average mass fractions in IBBN amount to as 
much as the solar system abundances. As see from Fig. [HI there are over-produced ele- 
ments around A = 150 (for rjf^^gf^ = 10"'') and A = 80 (for rjhigh = 10^^). It seems to be 
conflict with the chemical evolution of the universe. However, we show only the results of 
the upper-bounds on f^ — R diagram. Since /„ and R are free-parameters, over-production 
can be avoided by the adjustment of /„ and/or R. Figure |9] illustrates the mass fraction 
in r]f^^ , = 1.0 X 10"^ with various f^ — R sets. It is shown that the abundance pattern 
can be lower than the solar system abundance. Although we showed here only the result 
of ?7/jjg/i = 10~^ case, it is possible to avoid producing over-abundance in other parameters. 



Vhigh = 10 ^ and r^^^. ^ = 5x10 '^. If we put constraint on the fv — R plane from the heavy 



igh 

ill '^ o T-i /-| /n 

igh 

element observations, the limit of those parameters should be tightly. 

In our calculation, the radioactive nuclei are produced much in the high-density region. 
Especially, we should note that ^^Ni decays into ^^Fe (^®Ni— )■ ^®Co— >• ^^Fe), where the exis- 
tence of ^^Fe surely affects the process of the formation of the first generation stars. There- 
fore, it may be also necessary for IBBN to be constrained from the star formation scenarios, 
because opacity change due to IBBN will affect them. 

Recent observational signal of over-abundances of "^He mass fractions in globular clusters 
could motivate the IBBN scenario toward the detailed modeling. The over-abundances of 
^He are suggested to be in the range of 0.3 — 0.4 where estimated from the H-R diagram 



of the blue Main-Sequence of NGC2808 in Ref. 21j. If the origin of ^He in globular clusters 
is due to IBBN, rj must be greater than 10~^ in some regions during the epoch of BBN. 
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Then, the averaging procedure could be constrained from the more detailed observations 
of abundances. Since the history of changes in abundances has been investigated in detail 
through the chemical evolution of galaxies [35| , further plausible constrains on the averaging 
process should be studied in the next step. 

In our study, we ignore the diffusion effects. However, it is shown that the diffusion affects 



the primordial nucleosynthesis significantly 12| . Matsuura et al. [24| has estimated the size 



of the high-baryon density island to be 10 m - 10 m at the BBN epoch. The upper bound 

ower is from the analysis 



is obtained from the maximum angular resolution of CMB and the 

of comoving diffusion length of neutron and proton given in Ref . ll| . In our case, we can 

estimate the scale of the high-density and the effects of the diffusion from f^. 

The neutron diffusion effects can be discussed with use of the results obtained in the 
previous section by comparing the scale of the high-density region with the diffusion length. 
The present value of the Hubble length is H^^ = 1.28 x 10^^ m. We may estimate the scale of 
the high-density region from the Hubble length H~^ multiplied hy fv ■ From ranges of the 
volume fraction adopted in ^Vl 10~^ < fv < 10~^, we obtain the scale of the high-density 
regions at present epoch d^ as 2.7 x 10^'^ m. < d^ < 2.7 x 10^^ m. We can estimate the scale 
at redshift z from the relation d{z) = dQ{l + z)^^. As the result, we expect d at BBN era 
{z ~ 10^) as 2.7 X 10^^ m < d^^j^ < 2.7 x 10^^ m. We can say that the nucleon diffusion 
effects would be neglected because the diffusion length is much smaller than d. 

On the other hand, the high-density region is expected to be smaller than 10^^ m. It seems 
to be very bad that the upper bound of d is larger than the value as far as our two zone model 
is concerned. However, the high-density island cannot be observed directly, since we assume 
that the high- and low-density regions become homogeneous after the nucleosynthesis. 

Finally, distances between high density regions are difficult to derive without specific 
models beyond the two-zone model. We will plan to calculate the nucleosynthesis with the 
diffusion of abundances and/or more plausible averaging process included. 

Ackno^vledgments 

This work has been supported in part by a Grant-in- Aid for Scientific Research (18540279, 
19104006, 21540272) of the Ministry of Education, Culture, Sports, Science and Technology 
of Japan, and in part by a grant for Basic Science Research Projects from the Sumitomo 



13 



Foundation (No. 080933). 



[1] G. Steigman, Ann. Rev. Nucl. Part. Sci. 57, 463 (2007); 

F. locco, G. Mangano, G. Miele, O. Pisanti and P. D. Serpico, Phys. Rept. 472, 1 (2009) 
[2] V. Luridiana,A. Peimbert, M. Peimbert, k M. Cervino, Astrophys. J. 592, 846 (2003) 
[3] Y. L Izotov, T. X. Thuan and G. Stasinska, Astrophys. J. 662, 15 (2007) 
[4] Olive & Skillman, Astrophys. J., 617, 29-40 (2004) 
[5] D. Kirkman, D. Tytler, N. Suzuki, J. M. O'Meara and D. Lubin, Astrophys. J. Suppl. 149, 1 



(2003) [arXiv:astro-p h/0302006 1 . 
[6] M. Pettini, B. J. Zych, M. T. Murphy, A. Lewis, k C. C. Steidel, Mon. Not. R. Astron. Soc. 

391, 1499, (2008) 
[7] J. M. O'Meara, S. Buries, J. X. Prochaska, G. E. Prochter, R. A. Bernstein and K. M. Burgess, 

Astrophys. J. 649, L61 (2006) 
[8] C. Alcock, G.M. Fuller, and G.J. Mathews, Astrophys. J. 320, 439 (1987) 
[9] N. Terasawa and K. Sato, Phys. Rev. D 39, 2893 (1989) 



[10] K. Jedamzik, and J.B. Rehm, Phys. Rev. D64, 023510 (2001) |astro-ph/0101292 1 ; 

T. Rauscher, H. Applegate, J. Cowan, F. Thielmann, and M. Wiescher, Astrophys. J. 429, 

499 (1994). 
[11] J. H. Applegate, C. J. Hogan, and R. J. Scherrer, Phys. Rev. D35, 1151 (1987) 
[12] R. M. Malaney and W. A. Fowler, Astrophys. J 333, 14 (1988); 

J. H. Applegate, C. J. Hogan, R. J. Scherrer, Astrophys. J. 329, 572 (1988); 

N. Terasawa and K. Sato, Astrophys. J. 362, L.47 (1990); 

D. Thomas, D. N. Schramm, K.A. Olive, G. J. Mathews, B. S. Meyer, and B. D. Fields, 

Astrophys. J. 430, 291 (1994); 
[13] K. Jedamzik, G. M. Fuller, G. J. Mathews, and T. Kajino, Astrophys. J. 422, 423 (1994); 
[14] S. Matsuura, A. D. Dolgov, S. Nagataki and K. Sato, Prog. Theor. Phys. 112, 971 (2004) 
[15] G. M. Fuller, G. J. Mathews and C. R. Alcock, Phys. Rev. D 37, 1380 (1988); 

[16] H. Kurki-Suonio and R. A. Matzner, Phys.Rev. D39, 1046 (1989); 
H. Kurki-Suonio and R. A. Matzner, Phys.Rev. D42, 1047 (1990); 

14 



[17] C.L. Bennett, et al, Astrophys. J. Suppl. 148, 1 (2003) 

D. N. Spergel et al, Astrophys. J. Suppl. 170, 377 (2007) 

J. Dunkley et al. Astrophys. J. Suppl. 180, 306 (2009) 
[18] E. Komatsu et al. larXiv:1001.4538l [astro-ph.CO]. 
[19] Y. Juarez, R. Maiolino, R. Mujica, M. Pedani, S. Marinoni, T. Nagao, A. Marconi, &: E. Oliva, 

Astron. & Astrophys., 494, L25, (2009) 
[20] L. R. Bedin et al., Astrophys. J., 605, L125 (2004); 
[21] G. Piotto et al., Astrophys. J., 661 L53, (2007) 
[22] T. Rauscher, Phys. Rev. D 75, 068301 (2007) 
[23] S. Matsuura, S. I. Fujimoto, S. Nishimura, M. A. Hashimoto and K. Sato, Phys. Rev. D 72, 

123505 (2005) 
[24] S. Matsuura, S. I. Fujimoto, M. A. Hashimoto and K. Sato, Phys. Rev. D 75, 068302 (2007). 
[25] K. Jedamzik |astro-ph/9911242| . 



[27 
[28; 
[29; 
[30 
[31 

[32; 
[33 
[34; 
[35; 
[36 



[26] S. Fujimoto,M. Hashimoto, O. Koike,K. Aral, & R. Matsuba, Astrophys. J. 585, 418 (2003), 
O. Koike, M. Hashimoto, R. Kuromizu, & S. Fujimoto, Astrophys. J. 603, 592 (2004), 
S. Fujimoto, M. Hashimoto, K. Aral, & R. Matsuba, Astrophys. J. , 614, 847 (2004), 
S. Nishimura, K. Kotake, M. Hashimoto, S. Yamada, N. Nishimura, S. Fujimoto and K. Sato, 
Astrophys. J. 642, 410 (2006). 



M. Hashimoto & K. Aral, Physics Reports of Kumamoto University, 7, 47, (1985). 



B. Fields and S. Sarkar, arXiv:astro-ph/0601514. 
L. Kawano, FERMILAB-Pub-92/04-A 

E. Anders and N. Grevesse, Geochim. Cosmochim. Acta 53, 197 (1989). 

C. Angulo, M. Arnould, M. Rayet, P. Descouvemont, D. Baye, C. Leclercq-Willain, A. Coc, 
S. Barhoumi, P. Aguer, C. Rolfs, et al.. Nuclear Physics A 656, 3 (1999). 

M. Hashimoto, Progress of Theoretical Physics, 94, 663, (1995). 

K. Hagiwara et al [Particle Data Group], Phys. Rev. D 66, 010001 (2002). 

R. V. Wagoner, W. A. Fowler, & F. Hoyle, Astrophys. J. , 148, 3 (1967) 

M. E. Anderson, J. N. Bregman, S. C. Butler and C. R. Mullis, Astrophys. J. 698, 317 (2009) 

T. Moriya and T. Shigeyama, Phys. Rev. D 81, 043004 (2010) 



15 



o 
o 

CO 



CO 




J U U I U I I U l_ 



10 



10' 10" 10^ 10^ 
Time [sec] 



10^ 



(a) Vh 



igh 



1.02 X 10"^ 



o 

o 

03 



03 



10^ E 
10-2 - 

10-^ 
10-^ 

10"^ h 



10 
10 
10 



10 



-12 



14 




0° 10^ 10^ 10^ lO"^ 10^ 
Time [sec] 

(b) VMgh = 1-06 X 10-3 



Fig. 4: Time evolution of the mass fractions in high-density regions of (a) iJi^- , 
(b) VH^,^ = 1-06 X 10-3 



1.02 X 10"^ and 



high 



16 



1 



xlO^xlO^ 



temperature [K] 
1.0x10^ 4.0x10^ 



1.0x10 



10" 





10"" 






c^ 




-a 


10"^ 


3 




X) 


10"^ 


c^ 


n 




o 




^ 


10 


3 


10 


(i> 




n 






-12 




10 



10 



14 



10 




10 



10 10 

time [sec] 



10 



10 



Fig. 5: Time evolution of the neutron abundance in SBBN (r/ = 6.1 x 10 ^^) and IBBN in the 
high-density region. 



17 




C 11(29) 

[sli (2^) ^ 
Co(27) - 



Kr(3S) 

Br (35) ■- 

Se(34) - 

As|33) - 

Ge|32) p 

Ga(31)l 

Zn(30) — 



Fig. 6: Abundance distribution at T = 3.5 x 10^ K (t ~ 9 sec) in r]^-f^ = 5.3 x 10 ^. The gray 
regions are those of stable nuclei. 



18 




Ar(ia) 

CI (17) 
S (1^)-- 



Fig. 7: . Abundances distribution at T = 1.9 x 10^ K (t ~ 40 sec) in 77^. . = 5.3 x 10~^ 



10" 



10" 



10" 



.2 10 

CO Q 

^ 10"^ 

S -9 

fi 10 



10" 



10 



10 



11 



12 



w . 






ri,. , =1.0x10' 

-— x\,. =1.0x10'' 

ri,. =5.5x70" 

-■- ri,. =5.3x10"' 

— ■ ^,. =1.0x10"' 



1, , .*Ai ^ * 






100 150 

mass number 



200 



250 



Fig. 8: Comparison 
abundances 



of the averaged mass fractions in the two-zone model with the solar system 
30|] (indicated by dots). 



19 



Table. II: Mass fractions of heavy elements {A > 7) for three cases of ??;ij„/j — 10 ^,rihigh = 
5.33 X 10""^, and r/^.^^ ~ 10"^. 

/„ = 2.1 X 10-^ i? = 1.8 X 10'^ /„ = 7.0 X 10-^ i? = 9.3 X 10^ /^ = 5.7 x IQ-^ i? = 1.9 x 10^ 

{vMgK = 1-06 X 10-^) (r?,,^, = 5.33 x 10'^) {y.^^^ = 1.02 x IQ-^ ) 

element high average element high average element high average 

'Ni 1.25 X 10"^ 4.55 x 10"^ ^'^^Nd 1.99 x 10"^ 1.22 x 10"^ ^^^Nd 3.69 x lO"'^ 3.61 x 10"^ 
'Co 1.59 X 10"^ 5.79 x lO"'^ ^^Ni 1.37 x 10"^ 8.35 x lO"'^ ''°Ca 2.71 x lO"'^ 2.65 x 10"^ 
'Sr 1.06 X 10"^ 3.86 x lO"'^ ^''^Sm 1.03 x 10"^ 6.30 x lO"'^ ^^jyi^ 2.42 x lO"'^ 2.36 x 10"^ 



56 

57 
86 

^■^Sr 9.77 X 10-6 3.56 x lO"'^ ^^^Pm 8.91 x 10"^ 5.44 x 10""^ ^^^Eu 2.37 x lO"'^ 2.32 x 10" 
■^^Se 9.75 X 10-6 3.55 x 10"'^ ^''^Sm 8.25 x 10"^ 5.05 x lO"'^ ^^^Ce 1.93 x lO"'^ 1.89 x 10" 



'Se 2.93 X 10-6 1.07 x lO"'^ ^^^Pm 6.62 x 10"^ 4.05 x lO"'^ ^^Cr 1.55 x lO"'^ 1.51 x 10"^ 

^Sr 9.17 X 10-6 3.34 X lO"'^ ^''^Sm 5.24 x 10"^ 3. 20 x lO"'^ ^^^^e 1.11 x lO"'^ 1.09 x 10"^ 

-Kr 8.91 X 10-6 3.25 x 10"'^ ^^^p^^^ ^_i^ ^ ^Q-e 2.54 x lO"'^ ^6^^ ^ jq ^ iQ-^ 1.O8 x 10"^ 

Kr 7.80 X 10^6 2.34 x lO"'^ ^'''^Sm 3.99 x 10"^ 2.44 x lO"'^ ^^'^Nd 1.05 x lO"'^ 1.03 x 10"^ 

'Ge 7.67 X 10~6 2. 80 x lO^'^ ^^''Pm 3.66 x 10^6 2.24 x 10"^ 156e^ 9 44 ^ 10-8 9.22 x 10"^ 

'Kr 7.60 X 10-6 2.77 x lO"'^ ^^^p^^^ 345 ^ ;^q-6 2.12 x lO"'^ ^^^Nd 9.36 x 10"^ 9.15 x 10"^ 

'Kr 7.06 X 10-6 2.57 x lO"'^ ^'^^Nd 2.84 x 10"^ 1.74 x lO"'^ ^"^Fe 8.97 x 10"^ 8.77 x 10"^ 

'Kr 6.25 X 10^6 2.28 x lO"'^ ^''^Sm 2.67 x 10"^ 1.63 x lO"'^ ^^i^b 8.96 x 10"^ 8.52 x 10"^ 

'Ge 6.14 X 10^6 2.24 x lO"'^ ^''^Nd 2.25 x 10"^ 1.37 x lO"'^ ^^^La 8.80 x 10"^ 8.60 x 10"^ 

'Se 5.93 X 10^6 2. 16 x lO"'^ ^''^Sm 1.76 x 10"^ 1.07 x lO"'^ ^''N 8.74 x 10"^ 8.54 x 10"^ 

'Br 5.90 X 10^6 2.15 x lO"'^ ^^^Pm 1.16 x 10"^ 7^9 ^ lo-^ ^^Qr 8.56 x 10"^ 8.39 x 10"^ 

''Se 5.35 X 10-6 1.95 x lO"'^ ^^^Sm 9.88 x lO^'^ 6.04 x lO"'^ ^^^Ba 7.96 x 10"^ 7.77 x 10"^ 
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'Zr 4.41 X 10^6 iQi X 10-7 108^(1 592 x lO^^ 3.62 x lO"'^ ^62p,y g 34 ^ 10-8 6.68 x 10"^ 

'Rb 4.32 X 10^6 1 58 ^ 10"^ ^^^Eu 5.30 x lO^'^ 3.24 x lO^^^ ^^C 6.43 x 10"^ 6.28 x 10"^ 

'Rb 4.08 X 10^6 1 49 X 10-7 153eu 5.23 x lO^^^ 3.20 x 10"^ ^^O 6.30 x 10"^ 6.16 x 10"^ 

'Y 3.85 X 10-6 1.40 X lO^^^ ^^OCd 4.70 x lO^^^ 2.88 x 10"^ ^^^Gd 5.85 x 10"^ 5.71 x 10"^ 



89 

90 
85 
83 

88 ^ 

^^Zr 3.55 X 10^6 ^ 29 x lO"'^ ^^^Eu 3.79 x lO"'^ 2.32 x 10"^ ^^'^Cs 5.56 x 10"^ 5.43 x 10"^ 
73a« .-^ .^9 V 10-6 1.28 X lO-'^ i52j^^ 3.45x10^^7 2.11 



^As 3.52 X 10^6 1 28 x 10-'7 1^25^^ 345 ^ ^q^j 2.11 x 10^^ ^^■7Nd 3.96 x 10"^ 3.87 x 10"^° 

^^Ga 3.40 X 10-6 1.23 x 10-'7 ^'^^Ce 3.48 x 10-'7 2.12 x 10"^ ^^^Ro 3.77 x 10"^ 3.68 x 10"^ 

■^^Se 2.93 X 10^6 ^ Q7 ^ 10-'7 ^^Oj^^ 2.90 x lO"'^ 1.77 x 10"^ ^"^^Pr 3.11 x 10"^ 3.04 x 10"^ 

20 

^^Nb 2.90 X 10^6 1 06 x lO"'^ ^^^cj 2.89 x lO"'^ 1.77 x 10"^ ^^^Ce 3.00 x 10~^ 2.93 x 10"^° 

Y^ Y{ A\ Q in V in-4 1 iQ N/ in-5 Y^ vf a\ i 7s v in-4 7 on ^ in-6 Y^ y( A\ -i rk -^ in-6 q t;? ^, in-7 



o 



o 

CO 






10" 



10" 



10"° - 



10" 
10" 
10"" k 



10 



10 



10 



- + + 



•10 



•11 



-12 



+ 



y 



+4 



f„=1 .0x10"" R=1. 7x10^ 

L=1. 0x1 0;\R=1. 7x10^ 

fv=5.7x10"\ R=1.9x10^ 

Anders Grevesse1989 



\ 








100 150 

mass number 



200 



Fig. 9: Same as Fig. [U but i]^- , is fixed as 10 
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